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ABSTRACT 

The observation of the hot gas surrounding Sgr A* and a few other nearby galactic 
nuclei imply that electron and proton mean free paths are comparable to the gas 
capture radius. So, the hot accretion flows is likely to proceed under week-collision 
conditions. Hence, thermal conduction has been suggested as a possible mechanism 
by which the sufficient extra heating is provided in hot ADAF accretion disks. We 
consider the effects of thermal conduction in the presence of a toroidal magnetic field 
in an advection-dominated accretion flow around a compact object. For a steady-state 
structure of such accretion flows a set of self-similar solutions are presented. We find 
two types solutions which represent high and slow accretion rate. They have different 
behaviors with saturated thermal conduction parameter, (f>. 
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1 INTRODUCTION 

' An Advection-Dominated Accretion Flow (ADAF) is de- 
fined as one in which a large fraction of the viscously gener- 
ated heat is advected with the accreting gas, and only small 
fraction of the energy is radiated. Advection-dominated ac- 
cretion can occur in two different limits: 1- At very high 
mass accretion rates, radiation is trapped in the accreting 
gas because of the large optical depth and is advected with 
the flow. This limit of advection-dominated, typically occurs 
for mass accretion rates M > MEdd (the Eddington rate). 
2- At sufficiently low M, the accreting gas can become op- 
tically thin. The cooling time of the gas is then longer than 
the accreting time, and once again we have an ADAF. 

The thin advection-dominated accretion flow (ADAF) 
model has been investigated extensively since the end of 
1990s; (Ichimaru 1977, Narayan & Yi 1994, 1995a, 1995b; 
Abramowicz et al. 1995). This model has been used to in- 
terpret the spectra of black hole X-ray binaries in their qui- 
escent or low/hard state as an alternative to the Shapiro, 
Lightman & Eardly (1976, SLE) solutions. Since ADAFs 
have large radial velocities, and infalling matter carries the 
thermal energy into the black hole, advective energy trans- 
port can stabilize the thermal instability; thus ADAF models 
have been widely applied to explain observations of galactic 
black hole candidates (e.g., Narayan et al. 1996; Hameury et 
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al. 1997), the spectral transition of Cyg X-l (Esin et al. 1996) 
and multi-wavelength spectral properties of Sgr A (Narayan 
& Yi 1995b; Manmoto et al. 2000; Narayan et al. 1997). 
In addition, many ADAF-like models have been proposed 
including outflows, convection, etc. 

The theories for the structure and properties of hot 
accretion flows has been remained controversial. However, 
Narayan & Yi 1994 (hereafter NY1994) derived self-similar 
ADAF solutions which emphasized on the importance of the 
stabilizing role of the radial heat advection. Subsequent an- 
alytical work on the hot accretion flow has emphasized on 
the possibility of outflows, motivated by a positive Bernoulli 
constant (Narayan & Yi 1994, Fukue 2004). Also investi- 
gation was indicated a advection-dominated accretion flow, 
where the disk is convective in the radial direction (Narayan 
et al. 2000, Quatret & Gruzinov 2000). 

The diversity of the models indicate that modelling the 
hot accretion flows is a challenging and controversial prob- 
lem. We think that, one of the largely neglected physical 
phenomena is the thermal conduction while recent observa- 
tions of the hot accretion flow around active galactic nuclei 
indicated that it should be based on collision- less regime. 
So, thermal conduction probably has an important role in 
energy transport in the accreting materials in a hot accre- 
tion disk where they are nearly completely ionized. Recently, 
Medvedev & Narayan (2001) discovers a new type of accre- 
tion flow, a hot settling flow around a rapidly rotating neu- 
tron star. The flow is cooling-dominated and energetically 
similar to the SLE solution. The cooling dominated SLE so- 
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lution has been shown to be thermally unstable (Piran 1978, 
Wandel & Liang 1991; Narayan & Yi 1995a) and, hence, 
unlikely to exist in the nature. It has been shown that any 
accretion flow in which the heating balances cooling is ther- 
mally unstable, if the cooling is due to the bremsstrahlung 
emission (Shakura & Sanyev 1976; Piran 1978). But ADAF 
is known to be thermally stable (Narayan & Yi 1995a; Kato 
et al. 1996, 1997)therefore the cooling is week and the ther- 
mal energy of the flow is not radiated but advected with the 
gas. Since the SLE solution for non-ADAF disks is known 
to be thermally unstable, one might suspect that the new 
solution would also be unstable. However, due to the very 
high temperature of the accreting gas, thermal conduction 
is very strong and could suppress the thermal instability. 
So it should be important to consider the role of thermal 
conduction in a SLE solution. There is a branch of solu- 
tion for the problem of a smooth transition from an outer 
Shakura & Sanyev disk (SSD) to an inner ADAF disk around 
a compact object which include an additional mechanism of 
energy transport, thermal conduction. A number of authors 
showed that the SSD- ADAF transitions were realizable if an 
extra heat flux caused by thermal conduction was invoked ei- 
ther in the radial direction (Honma 1996; Manmoto & Kato 
2000; Gracia et al. 2003) or in the vertical direction (Meyer 
& Meyer-Hofmeister 1994; Meyer, Liu, & Meyer-Hofmeister 
2000). 

The weakly-collisional nature of ADAFs has been noted 
previously (Mahadevan & Quataret 1997). But a few au- 
thors tried to study the role of turbulent heat transport in 
ADAF disks ( Honma 1996, Manmoto et al. 2000). Since 
the thermal conduction acts to oppose the formation of the 
temperature gradient that causes it, that might expect that 
the temperature and density profile in a thermal conduct- 
ing disk should be different from the case in which the disk 
is not under the influence of thermal conduction. Recently, 
Tanaka & Menou (2006), studied the effect of saturated ther- 
mal conduction on optically thin ADAFs using an extension 
of self-similar solution of NY1994. Their solutions suggest 
that the thermal conduction may have an important role in 
the dynamical behavior of the hot accretion flows and prob- 
ably is an important factor to understand the physics of hot 
accretion disks. 

The aim of this work is to consider the possibility of 
the thermal conduction in the presence of toroidal magnetic 
field, which has been largely neglected ingredient before, 
could affect the global properties of the hot accretion flows 
substantially. The tangled magnetic fields in accretion flows 
would likely reduce the effective mean free paths of particles. 
The magnitude of this reduction, which will depends on the 
field geometry, is still unknown. In this manuscript we will 
investigate the effect of thermal conduction on the physical 
structure of ADAF-like accretion flow around a black hole 
in the presence of a toroidal magnetic field. 



2 THE BASIC EQUATIONS 

Let us consider a gaseous disk rotating around a 
Schwarzschild black hole of mass M. The disk is assumed to 
be in an advection-dominated regime. We assume a steady 
axi-symmetry accretion flow ( = ^ = 0) and a geometri- 
cally thin disk. In cylindrical coordinates (r, cj>, z), we ver- 



tically integrate the flow equations. Also, we suppose that 
all flow variables are only a function of r. We ignored the 
relativistic effect and we use Newtonian gravity. The disk is 
supposed to be turbulent and possesses an effective turbu- 
lent viscous. We adopt a-prescription for viscosity of rotat- 
ing gas. We have assumed that the generated energy due to 
viscosity dissipation and heat conduction into the volume 
are balanced by advection cooling. The magnetic field was 
considered with toroidal configurations. 

We can describe the accretion flows by the fundamental 
governing equations which are written by the equation of 
continuity which is integrated in the vertical direction is 
expressed for presence purpose as 

1 -f r (rXv r )=2pH, (1) 

where v r is the accretion velocity, p the mass-loss rate per 
unit volume, H the disk half-thickness, and E the surface 
density, which is defined as E = 2pH. 

The equation of motion in the radial direction is: 

dv r V% GM Id 2 C \ 1 d 2 

^ = -~ — -E* (SCs) -T-2E^ (Sc ^' (2) 
where v^, c s , ca, are the rotation velocity of gas disk, sound 
speed and Alfven speed, respectively. Sound speed is defied 
as c 2 s — , pgas being the gas pressure and Alfven veloc- 
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ity is defined as c\ — — Pm p ag , where p mag being the 
magnetic pressure. The first term in the right hand side of 
this equation implies the centrifugal force while the third 
term implies the pressure gradient force. The fourth and 
fifth terms represent the magnetic force. 

By integrating the z-component of the momentum equa- 
tion along the z-axis, we have hydrostatic balance in the 
vertical direction as: 

^ 2 ^[1 + ^) 2 ] = (1 + /3) C *, (3) 

where f3 = p ™ ag — I(£A) 2 which indicates the importance 
of magnetic field pressure compared to gas pressure. We as- 
sume that this value is constant through the disk, but we 
will show the dynamical properties of the disk for different 
values of f3. 

For writing the angular momentum equations we need 
to choose proper mechanism for viscosity. We adopt a- 
prescription for viscosity. The r<^-component of the viscosity 
stress tensor is proportional to pressure (Shakura & Sanyev 
1973). 

dQ 

ir4> = vr — = a P> 
dr 

where rj = pv is the viscosity, v is the kinetic viscosity and 
a is a viscous parameter. Here we can see that the viscosity 
stress tensor is proportional to the gas pressure or to the 
total (gas and magnetic) pressure. So, we can choose two 
cases: Case 1: when the pressure is assumed to be the gas 
pressure (thermal pressure). Case 2: when the pressure is 
assumed to be the magnetic pressure plus the gas pressure. 
For simplicity in this investigation we use the case 1, we 
adopt the form: 

v = fi^'ttf?) = aC 3 H, (4) 
P 

This form is so-called a-prcscription where a is a constant 
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less than unity (Shukura & Sanyev 1973). So, the angular 
momentum transfer equation is: 



_ d , . d ,aEc?r 3 dQ. 
rZv r -(rv,) = -(-^—- J 



(5) 



is Keplerian 



where Q is angular velocity and Q, k = 
angular velocity. 

Now we can write the energy equation considering cool- 
ing and heating processes in an ADAF. The energy equation 
is: 

Eiv dc? Sc^ d 



7" 



1 dr 



-{vVr) — Qvis Qrad H~ Qcond-, 



(6) 



where Q„i 3 , Qcond and Q ra d are the locally released ener- 
gies due to viscous dissipation, transported by the thermal 
conduction and radiative cooling rate, respectively. 

The viscosity generates thermal energy through the dif- 
ferential rotation of turbulently moving gas, and the viscous 
heating rate is expressed by: 



Qvi 



Q 



■ rT r 



dtt 

' dr 



ac 2 B r 2 , dn , 
Slk ( dr' 



The radiated cooling rate is expressed as: 
8acT 4 



2>npH 

where k is the opacity of the rotating gas. We assumed that 
the disk to be radiation-pressure-dominated with opacity 
due to the electron scattering only. In the right hand side of 
the energy equation we have: 

Qvis Qrad "t - Qcond — Qadv 

where Q a dv represents the advective transport of energy 
and is defined as the difference between the viscous heating 
rate, Qvis, and radiative cooling rate, Q ra d, plus the energy 
transported by conduction, Qcond- We employ the parame- 
ter /= 1 — ® ) rad to measure the high degree to energy arc 
transported by radiation. When / ~ 1 radiation can be neg- 
ligible so it is in advection dominated state while for small 
/ disk is in radiation dominated case. So we can rearrange 
the right hand side of Eq.6 to fQ V is + Qcond, where / ^ 1. 
In general, / varies with r and depends on the details of 
heating and cooling processes. For simplicity, it is assumed 
a constant. 

The energy transport by thermal conduction, Q CO nd, is: 



Qc 



2H d I V X 

3-( rF ")' 

r dr 



(7) 



adopted from the formulation of Cowie & McKee (1977) and 
Tanaka & Menou (2006). The saturated conduction flux is 
then written as F s = 5<j)p(? a , where (f> is the saturated con- 
stant (presumably < 1), p is the gas mass density and c s is its 
sound speed. In order to implement the thermal conductiv- 
ity as correctly as possible it is essential to know whether the 
mean free paths is less than (or comparable to) the length 
scale of the temperature gradient. For electron mean free 
paths which are greater than the length scale of temperature 
gradient, the thermal conductivity is said to be 'saturated'. 
Because it is a saturated flux, it no longer explicitly depends 
on the magnitude of the temperature gradient but only does 
on the direction of this gradient. The heat will flow outward 
in a hot accretion flow with a near-virial temperature profile, 
hence the positive sign was adopted for F s . 



Finally, since we consider the toroidal component for 
the global magnetic field of central stars, the induction equa- 
tion with field scape can be written as: 



(8) 



where is the field escaping or crating rate due to the 
magnetic instability or dynamo effect. We can rewrite this 
equation as: 



dr 
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(9) 



Now we have a set of MHD equations which describe 
the dynamical behavior of ADAF flows. The solution of these 
equations gives us the dynamical behavior of the disk, which 
strongly depends on the viscosity, magnetic field strength, 
thermal conduction and advection rate of energy transport. 



3 SELF-SIMILAR SOLUTIONS FOR ADAFS 
WITH SATURATED CONDUCTION 

The self-similar solution can not be able to describe the 
global behavior of the accretion flow, because in this method 
there are not boundary conditions which have been taken 
into account. However as long as we are not interested in 
the solutions near the boundaries, such solutions describe 
correctly the true and useful solutions asymptotically at in- 
termediate areas. 

Narayan & Yi (1994) simplified 2D axi-symmetric 
steady state ADAF problem by assuming that the dynami- 
cal variables of the flow have a power law dependence on r. 
This allowed them to evaluate directly all radial derivations 
in the equations. Further, they eliminated 6 by considering 
a height-integrated set of equations. The equation has an 
analytical self-similar solution which depends only on three 
parameters: the viscous parameter a, the ratio of specific 
heats of the accreting gas 7, /. In our case, we have a mag- 
netic field in the structure of the ADAF. So, we have an 
extra free parameter /3, which introduces the magnetic field 
strength in our self-similar solution. 

We assume that each physical quantity can be expressed 
as a power law of the radial distance, r" , where power index 
v is determined for each physical quantities self-consistently. 
The solutions are: 



v r = — Cia 



GM 



C 2 
1 V 



GM 
r 

GM 



c s = - = Cz — 
p r 



S = Sor _ 5, 
p = p r' 3 , 
Bs = B r~^, 



GM 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 
(16) 
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where CI, C2 and C3 are coefficients that we will determine 
later, po and Bo are constants, which provide convenient 
units with which the equations can be written in the non- 
dimensional form. 

Using these solutions from the continuity, the momen- 
tum, angular momentum, hydrostatic, energy and induction 
equations, we can obtain the following system of dimension- 
less equations, solved for CI, C2, C3, / and cS: 



p = 0, 

ia 2 C 2 + C 2 2 -l + i[3-/3]C 3 = 0, 



Ci 

H 

r 

cl 

B 



3C 3 
2 : 



= V(l + /9)Cs, 



3-7 2 40 -1 

7-19/ 9a/, ft 

•' V 2 



5CiaGM 



47tB c 



/3C 3 



V(1+/3)C 3 



(17) 
(18) 

(19) 

(20) 

(21) 

(22) 



This solution tends to the solution presented by Akizuki & 
Fukue (2006) for = 0. Our solution is compatible with 
the standard solution of ADAF, Narayan & Yi 1994, hav- 
ing no mass loss (|17|l . Actually, radial dependence of surface 
density, (| 14f> . is a free parameter. For finding a physical so- 
lution, it should be more than -1. We chose it equal to 
to have a physical solution which is compatible with the 
standard solution in the case there is no any wind or mass 
loss effects. When it is more than — i we have mass loss 
in the disks. Fukue (2004) solved the standard disk for the 
case that radial dependence of surface density is equal to 2. 
He showed that the global behavior of the mass loss disk is 
quite similar to advective disks. In additions, we can easily 
see from (|20|l that the disk thickness becomes large due to 
magnetic pressure. 

In the case of very small a, we have Bo = which means 
that the creation or escape, so the toroidal component of 
magnetic field are balanced each other. In this case, the disk 
supports itself with rotation, gas and magnetic pressure. In 
the case of a finite a, the effect of injection of the toroidal 
component of magnetic field, can moderate the dynamical 
behavior of the disks. 

The above equations describe self-similar behavior of 
the optically thick advection dominated accretion disk with 
saturated thermal conduction in the presence of toroidal 
magnetic field. In a finite a, we will solve these equations. 
After some algebraic manipulations we can find a fourth or- 
der equation for Ci as follow: 

D 2 C\ + 2BDC\ + (B 2 - 2D)Cl - {A 2 + 2B)d + 1 = 0, (23) 
where 

D = -a 2 , 
2 ' 
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Figure 1. Profile of Ci; coefficient of infall velocity as a function 
of / ( which represent the importance of radiation trans- 
port). This profile illustrates the solutions of forth order equation 
for Ci in a given range of parameters. The solid line represents 
the low accretion rate solution while the dashed line represents 
the high accretion rate solutions. 



This Algebraic equation shows that the variable Ci 
which determines the behavior of radial velocity depends 
only on the a, <f>, (5 and /. Other flows quantity such as C2 
and C3 can be obtained easily from Ci. 

This equation is plotted in Fig.l for the given parame- 
ters. We have found that it has two real roots for this range 
of parameter space with two distinct behaviors. One of them 
represents low accretion rate and other one represent high 
accretion rate where both of them are sub-Keplerian. The 
parameters of the model are the ratio of specific heats 7, 
the standard viscous parameter a, the radiation transport 
parameter / and the degree of magnetic pressure to gas pres- 
sure p. 

Figure 1 shows the coefficient Ci, which represents the 
behavior of the radial flows of accreting materials. Although 
in ADAF model the radial inf all ve locity is generally slower 
than the Keplerian speed {~\/ ^-), it becomes large with /. 
As the level of advection is increased (by increasing the /), 
the accreting materials increase their inflow speed in both 
two distinct solutions. This variation is around 2 or 3 times 
for different range of /. This dependence is consistent with 
the usual ADAF solutions (Akizuki & Fukue 2006, Tanaka 
& Menou 2006). 

In Fig. 2, we investigate the role of saturated thermal 
conduction on the radial flow for the two above solutions. 
These solutions have different behaviors as we increase the 
role of thermal conduction by adding <f). The results for 
<f> = 0.001 to 0.006 are shown in this manuscript. Tanaka 
& Meneu (2006) have shown that for a very small (j> their 
solutions coincide the original ID ADAF solutions; but by 
adding saturated conduction parameter, <fi, the effect of ther- 
mal conduction can be better seen when we approach to 
~ 0.001 — 0.01. So we have plot our solution in this range. 
As it can bee seen in Figure 2 (upper) in the low accretion 
rate solution, when we increase <f> the radial flow decreased. 
But for the high accretion rate solution the radial flow is 
increases by increasing the <f> (Fig. 2 lower). Increasing the 
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havior is inversed. The behavior of high accretion solution is 
similar to the solution presented by Ghanbari et al. (2007). 
In this case, Increasing the viscous parameter is correspond- 
ing to the increase of the heating mechanism, so in a fixed 
advection regime, there is more energy to advect into the 
central star. 

Fig. 4 shows the coefficient C\ as a function of / for dif- 
ferent values of (3. By adding (3, which indicates the role of 
magnetic field in the dynamics of accretion disks, we will see 
that the radial flow increases in both solutions. On the other 
hand, the radial infall velocity increases when the toroidal 
magnetic field becomes large. This is due to the magnetic 
tension terms, which dominates the magnetic pressure term 
in the radial momentum equation that assists the radial in- 
fall motion. 
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Figure 2. Numerical coefficient Ci as a function / for different 
values of thermal conduction parameter, <j>. (Up) for low accretion 
rate solution and (Down) for high accretion rate solution. 

saturated conduction in high accretion solution will change 
the behavior of the radial flow in different range of advec- 
tion parameter. The effect of thermal conduction for low / 
is larger than compare to / ~ 1. The thermal conduction in 
low / causes a more radial flow. In the high accretion solu- 
tion, by adding the role of advection (adding /), the effect 
of thermal conduction was reduced. But, for high accretion 
solutions, the global effect of thermal conduction increases 
the radial flow, whereas in a low accretion solution, this role 
is quite suppressed. 

To Show the behavior of the solutions, self-similar ra- 
dial velocity is demonstrated in Fig. 3 for different values / 
with variation of viscosity parameters. The most accepted 
value for a is less than 0.1. There are probably, more signif- 
icantly variations when a exceeds above 0.1. Recently, King 
et al. (2007) assert that in a thin and fully ionized disk, 
the best observational evidence suggests a typical range of 
a ~ 0.1 — 0.4 where relevant numerical simulations tend to 
drive the estimates for a which are one order of magnitude 
smaller. However, such large values of a are unlikely (eg. 
Narayan, Loab, & Kumar 1994 , Hawley et al. 1994), so we 
have not explored this region a. Here we investigate these 
solutions with different values of viscosity less than 0.1. In 
the low accretion solution (Fig. 3 upper), when we increase 
the viscosity parameter, we see radial inflow increases. But 
in the case of high accretion solution (Fig. 3 lower) this be- 



4 CONCLUSION 

In this paper we have studied an accretion disk around a 
black hole in an advection-dominated regime in the presence 
of toroidal magnetic field with thermal conduction. We have 
presented the results of self-similar solutions to show the ef- 
fect of the viscosity and thermal conduction on magnetically 
driven accretion flows from a flow threaded by toroidal mag- 
netic fields. The only serious approximation we have here 
made is the use of an isotropic a viscosity. Attention has 
restricted to flow accretions in which self-gravitation is neg- 
ligible. Considering weekly collisional nature of hot accretion 
flows, a saturated form of thermal conduction has been used 
as a possible mechanism for transporting the energy was pro- 
duced by viscous heating. In a sense, we have accounted for 
this possibility by allowing the saturated conduction con- 
stant, <f), to vary in our solutions. 

We have found a self-similar solution using similarity 
technique in analogy to self-similar solutions by Narayan 
& Yi (1994) and Akizuki & Fukue (2006). The self-similar 
solution is valid in the intermediate region of the advection- 
dominated disks. The toroidal geometry for magnetic field is 
a good approximation for magnetic structure in this region. 
The disk structure, is characterized by the coefficient Ci, 
slightly depends on the strength of magnetic field. The ge- 
ometrical thickness of the disks depends the magnetic field 
strength; while the structure of accretion flow will be modi- 
fied strongly by viscous parameter a and saturated thermal 
conduction parameter, (f>. 

We have found two solutions with different infall veloc- 
ity, where they have different behavior with thermal con- 
duction. In the high accretion rate solution, when the level 
of saturated conduction increases, the radial infall velocity 
also increases, while for the low accretion rate, the solution 
is quite different. 

The presence of magnetic field with toroidal geometry 
will affect the role of thermal conduction. In both solutions, 
we have found that a large magnetic field causes a more 
radial infall velocity. 

There are some limitations in our solutions. One of them 
is that the self-similar hot accretion flow with conduction is 
1-temperature structure. If we use a 2-temperatures struc- 
ture for the ions and electrons in the disks, it is expected 
that the ions and electrons temperatures will decouple in 
the inner regions, which will modify the role of conduction. 
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Figure 3. Numerical coefficient C'i as a function / for differ- 
ent values of viscous parameter, a. (Up) for low accretion rate 
solution and (Down) for high accretion rate solution. 

The other limitation of our solution is the anisotropic char- 
acter of conduction in the presence of magnetic field. Balbus 
(2001) has argued that the dynamical structure of the hot 
flows could be affected by the anisotropic character of ther- 
mal conduction in the presence of magnetic field. 

However, our results clearly improve the physics of 
advection-dominated accretion flow around a black hole. It 
is important to investigate the effect of thermal conduction 
on the physical structure of hot flow around a black hole. 
We developed Narayan & Yi (1994) standard solutions to a 
more realistic model of ADAFs by adding the magnetic field 
on the structure of disks and thermal conduction as a mech- 
anism to transport energy in radial direction. In the future 
studies, we plan to improve our model with a more realistic 
viscous model, Case 2, and to follow the effect of thermal 
conduction to other physical parameters of the disks. Also, 
we are going to develop our model in a 2-D self-similar so- 
lution for ADAF around a black hole. Several developments 
can be investigated to reach a much more realistic descrip- 
tion for the physics of hot accretion disks around a magne- 
tized compact object. 
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A simple Science Template 



John Smith, 1 * Jane Doe, 1 Joe Scientist 2 

department of Chemistry, University of Wherever, 
An Unknown Address, Wherever, ST 00000, USA 
2 Another Unknown Address, Palookaville, ST 99999, USA 

To whom correspondence should be addressed; E-mail: jsmith@wherever.edu. 

This document presents a number of hints about how to set up your Science 
paper in WTjjX . We provide a template file, scif ile . tex, that you can use 
to set up the ETgX source for your article. An example of the style is the special 
{sciabstract} environment used to set up the abstract you see here. 

Introduction 

In this file, we present some tips and sample mark-up to assure your ETgX file of the smoothest 
possible journey from review manuscript to published Science paper. We focus here particularly 
on issues related to style files, citation, and math, tables, and figures, as those tend to be the 
biggest sticking points. Please use the source file for this document, scifile.tex, as a 
template for your manuscript, cutting and pasting your content into the file at the appropriate 
places. 

Science's publication workflow relies on Microsoft Word97. To translate ETgX files into 
Word97, we use an intermediate MS-DOS routine (?) that converts the T^X source into HTML. 
The routine is generally robust, but it works best if the source document is clean ETgX without 
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a significant freight of local macros or . sty files. Use of the source file scif ile . tex as a 
template, and calling only the .sty and . bst files specifically mentioned here, will generate 
a manuscript that should be eminently reviewable, and yet will allow your paper to proceed 
quickly into our production flow upon acceptance (?). 

Formatting Citations 

Citations can be handled in one of three ways. The most straightforward (albeit labor-intensive) 
would be to hardwire your citations into your ETgX source, as you would if you were using an 
ordinary word processor. Thus, your code might look something like this: 

However, this record of the solar nebula may have been 
partly erased by the complex history of the meteorite 
parent bodies, which includes collision-induced shock, 
thermal metamorphism, and aqueous alteration 
({\it 1, 2, 5 — 7\/}) . 

Compiled, the last two lines of the code above, of course, would give notecalls in Science style: 

. . . thermal metamorphism, and aqueous alteration (1, 2, 5-7). 

Under the same logic, the author could set up his or her reference list as a simple enumera- 
tion, 

{ \bf References and Notes} 
\begin { enumerate } 

\item G. Gamow, {\it The Constitution of Atomic Nuclei 
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